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Abstract 

We identify scenarios not ruled out by LEP data in which NMSSM Higgs detection 
at the LHC will be particularly challenging. We first review the ‘no-lose’ theorem for 
Higgs discovery at the LHC that applies if Higgs bosons do not decay to other Higgs 
bosons — namely, with L = 300 fb~^, there is always one or more ‘standard’ Higgs 
detection channel with at least a 5a signal. However, we provide examples of no- 
Higgs-to-Higgs cases for which all the standard signals are no larger than 7a implying 
that if the available L is smaller or the simulations performed by ATLAS and CMS 
turn out to be overly optimistic, all standard Higgs signals could fall below 5a even in 
the no-Higgs-to-Higgs part of NMSSM parameter space. In the vast bulk of NMSSM 
parameter space, there will be Higgs-to-Higgs decays. We show that when such decays 
are present it is possible for all the standard detection channels to have very small 
significance. In most such cases, the only strongly produced Higgs boson is one with 
fairly SM-like couplings that decays to two lighter Higgs bosons (either a pair of the 
lightest CP-even Higgs bosons, or, in the largest part of parameter space, a pair of 
the lightest CP-odd Higgs bosons). A number of representative bench-mark scenarios 
of this type are delineated in detail and implications for Higgs discovery at various 
colliders are discussed. 
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1 Introduction 


One of the most attractive supersymmetric models is the Next to Minimal Supersymmetric 
Standard Model (NMSSM) [1] which extends the MSSM by the introduction of just one sin¬ 
glet superheld, S. When the scalar component of S acquires a TeV scale vacuum expectation 
value (a very natural result in the context of the model), the superpotential term \SHuHd 
generates an effective fiHuHd interaction for the Higgs doublet superhelds with /i = A (S'). 
Such a term is essential for acceptable phenomenology. No other SUSY model generates 
this crucial component of the superpotential in as natural a fashion. We also note that the 
LEP limits on Higgs bosons imply that the MSSM must be very highly hne-tuned, whereas 
in the NMSSM parameter choices consistent with LEP limits can be found that have very 
low hue-tuning [2,3]. Thus, the phenomenological implications of the NMSSM at future 
accelerators should be considered very seriously. 

In the NMSSM, the h, H, A, Higgs sector of the MSSM is extended so that there are 
three CP-even Higgs bosons {hi^ 2 , 3 , ^hi < ^h 2 < ^hs), two CP-odd Higgs bosons (ai^ 2 , 
niai < maa) (we assume that CP is not violated in the Higgs sector) and a charged Higgs 
pair (h^). Hence, the Higgs phenomenology in the NMSSM can differ signihcantly from the 
one in the MSSM (see refs. [4-9] for recent studies). 

Our focus will be on NMSSM Higgs discovery at the LHC. An important question is then 
the extent to which the no-lose theorem for MSSM Higgs boson discovery at the LHC (see 
refs. [10,11] for CMS and ATLAS plots, respectively) is retained when going to the NMSSM; 
i.e. is the LHC guaranteed to hnd at least one of the hi, 2 , 3 , oig, 

We will hnd that it is not currently possible to claim a no-lose theorem for Higgs discov¬ 
ery in the NMSSM. This is due to the importance of Higgs-to-Higgs decays in the NMSSM. 
Indeed, the no-lose theorem for MSSM Higgs boson discovery at the LHC is based on Higgs 
decay modes (hereafter referred to as ‘standard’ modes) other than Higgs-to-Higgs decays.^ 
The importance of such decays was hrst noted in [12] and later pursued in [13,14]. Corre¬ 
spondingly, the parameter space of the NMSSM can be decomposed into the following three 
regions: 

a) An (actually fairly small) region where, for kinematical reasons, Higgs-to-Higgs decays 
are forbidden. Here, Higgs detection in the NMSSM proceeds via the standard discovery 
modes, with possibly reduced couplings and altered branching ratios with respect to the 
MSSM. In a hrst exploration of this part of the NMSSM parameter space [12], signihcant 
regions were found such that the LHC would not detect any of the NMSSM Higgs bosons. 
Since then, however, there have been improvements in many of the detection modes (and 
the addition of new ones). As a result [4], if the neutral NMSSM Higgs bosons do not decay 
to other Higgs bosons, then the LHC is guaranteed to discover at least one of them for an 
integrated luminosity of L = 300 fb“^ at both the ATLAS and the CMS detectors. 

b) The largest region of the NMSSM parameter space is the part where Higgs-to-Higgs 
decays are kinematically possible, but where the standard discovery modes are still sufficient 

^Higgs-to-Higgs decays do not create a problem for the CP-conserving MSSM no-lose theorem due to the 
constrained nature of the MSSM Higgs sector. Relations among the MSSM Higgs boson masses are such 
that Higgs pair decays are only possible if ruA is quite small. In this part of parameter space, the H is 
SM-like and H — > AA decays can be dominant. However, when ttia is small, the h also has small mass and 
the Z hA coupling is large. As a result, Z hA pair production would have been detected at LEP. 


2 



for the detection of at least one Higgs boson at the LHC. 

c) For a small part of the NMSSM parameter space Higgs-to-Higgs decays are dominant 
for the Higgs bosons with snbstantial prodnction cross sections, and the standard discovery 
modes do not yield a 5cr signal (even for integrated Inminosity of L = 300 fb“^ and after 
combining modes) for any of the Higgs bosons. In Refs. [6,7], we presented a selection of 
benchmark points with these characteristics. However, since then the expected statistical 
signihcances at the LHC in the standard discovery channels have improved and some of these 
points wonld now give a 5a signal. (On the other hand, in [6,7] we were somewhat optimistic 
with respect to LEP limits on Higgs bosons with masses below 115 GeV and nnconventional 
decay modes.) 

In this paper we repeat these stndies, npdating the LEP constraints and the expected 
statistical signihcances for the standard discovery modes at the LHC. Once again, we hnd 
a region in the NMSSM parameter space of type c) above. In section EJ we present new 
benchmark points for which the primary decaying nentral Higgs boson {hn) has strong 
conpling to gange bosons and has mass in the range [90 GeV, 150 GeV] bnt decays almost 
entirely to a pair of even lighter secondary Higgs states (h^hi). Both the primary and 
secondary Higgs bosons will have escaped LEP searches and will be impossible to observe at 
the LHG in the standard modes. The benchmark points presented are chosen to represent a 
range of and uih^ possibilities and a variety of possible decays. 

The ontline of the paper is as follows: In section |21 in preparation for onr discussions, we 
dehne the NMSSM model and its parameters. There, we also review the program NMHDE- 
GAY [9] employed for this study, and specify our precise scanning procedures for the Higgs 
discovery studies. 

In sectional we review the conclusions of [4] regarding the above region a) of the NMSSM 
parameter space. These remain unchanged: assuming an integrated luminosity of L = 
300 fb“^, one can establish a no-lose theorem for the very restricted part of parameter space 
where there are no decays of neutral Higgs bosons to other Higgs bosons. The statistical 
signihcances as a function of the charged Higgs mass, and the properties of two relatively 
difficult points in this region of parameter space (but still with a 5cr signal in at least one of 
the standard discovery modes) are presented. 

In section 15 we discuss general properties of points for which one or more Higgs-to-Higgs 
decays are allowed and, as a result, discovery of a Higgs boson in one of the standard modes 
is not possible. We present eight new benchmark points, discuss their properties, and show 
how Higgs-to-Higgs decays can lead to very small signals in all the usual LHG Higgs discovery 
channels for these points. 

In section |5 we will discuss the nature and detectability of the collider signals for the 
Higgs pair decay modes, especially focusing on the difficulties at hadron colliders such as 
the Tevatron and LHG. Notably we propose that the LHG may be able to detect Higgs-pair 
hnal states using the WW hu ^ h^hi —>• jjT^T~ production/decay mode, and discuss 
its properties and possible cuts. Final conclusions are given in section |5 
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2 The model, the NMHDECAY program and the scan¬ 
ning procedures 

We consider the simplest version of the NMSSM [1], where the term HH 1 H 2 in the superpo¬ 
tential of the MSSM is replaced by (we use the notation A for the superheld and A for its 
scalar component held) 

XH.H^S + ^ ( 1 ) 

so that the superpotential is scale invariant. The associated trilinear soft terms are 

XA.SH^H, + '^A^S^. ( 2 ) 

The hnal two Higgs-sector input parameters are 

tan/3= = X (S) . (3) 

These, along with mz, can be viewed as determining the three Susy breaking masses squared 
for Hui Hd and S appearing in the soft-SuSY-breaking terms 

ml.Hf + ( 4 ) 

through the three minimization equations of the scalar potential. Thus, we make no assump¬ 
tion of “universal” soft terms. 

In short, as compared two independent parameters in the Higgs sector of the MSSM 
(often chosen as tan/? and M^), the Higgs sector of the NMSSM is described by the six 
parameters 

X , K , Ax , A^, tan/3 , /iefr • (5) 

We will choose sign conventions for the helds such that A and tan/3 are positive, while k, 
Ax, Hk and Hes should be allowed to have either sign. We will perform a scan over these 
parameters using the publicly available program NMHDECAY [9]. For any choice of the 
above parameters and other soft-SuSY-breaking parameters that affect radiative corrections 
and Higgs decays, NMHDECAY performs the following tasks: 

1. It computes the masses and couplings of all the physical Higgs and sparticle states. We 
only retain points for which all Higgs and squark/slepton masses-squared are positive. 

2. It checks whether the running Yukawa couplings encounter a Landau singularity below 
the GUT scale. In our scans, we eliminate such cases. 

3. NMHDECAY checks whether the physical minimum (with all vevs non-zero) of the 
scalar potential is deeper than the local unphysical minima with vanishing (Hu), (Hd) 
or (S). We keep only parameter choices for which the minimum with all vevs non-zero 
is the true minimum. 

4 . It computes the branching ratios into two particle hnal states (including charginos and 
neutralinos) of all Higgs particles. Currently, squark and slepton decays of the Higgs 
are not computed. 
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5. It checks whether the Higgs masses and couplings violate any bounds from negative 
Higgs searches at LEP, including many quite unconventional channels that are relevant 
for the NMSSM Higgs sector. It also checks the bound on the invisible Z width 
(possibly violated for light neutralinos). In addition, NMHDECAY checks the bounds 
on the lightest chargino and on neutralino pair production. Parameter choices that 
conflict with LEP bounds are eliminated in the scans discussed below. 

Points that pass the requirements of items 1 through 3 above dehne the set of “physically 
acceptable” parameter choices. Our scans will be for randomly chosen parameter values in 
the following ranges: 

10-^ < A < 0.75; -0.65 <k< 0.65; 1.6 < tan/3 < 54 

-1 TeV < /ieff, Ha, < +1 TeV. ( 6 ) 


In the gaugino sector, we chose M 2 = 1 TeV (at low scales). We assume universal 
gaugino masses at the coupling constant unihcation scale, leading to Mi ~ 500 GeV and 
M 3 3 TeV. Thus, the lightest neutralino can only be signihcantly lighter than 500 GeV if 
it is mainly singlino or (when yUgfr is relatively small) higgsino. For the chosen Mi 2,3 values, 
LHG detection of the gauginos will be quite difficult and decay of Higgs bosons to gauginos, 
including the invisible XiXi channel, will in most cases be negligible. 

Gurrent lower limits from LEP and the Tevatron imply that squarks and sleptons must be 
at least moderately heavy. As a result, the Higgs bosons with substantial WWjZZ coupling 
(which are predicted to have masses below about 150 GeV) cannot decay to squarks and 
sleptons. We will choose squark/slepton parameters in the TeV range. In this case decays to 
squarks and sleptons are unimportant or absent also for the heavier Higgs states. Specihcally, 
we choose mq = mu = mo = rnu = mE = Msusy = 1 TeV for the soft-SuSY-breaking 
masses for all generations. This means that squarks and sleptons will be at the edge of the 
LHG discovery reach and that Higgs boson detection might be the only new physics signal.' 

For the trilinear soft-SuSY-breaking squark parameters, we choose Au = A^ = 1-5 TeV 
for all generations, including the third generation. We recall that the light Higgs mass is 
maximized by choosing the parameter At so that 


Xt = 


At — fteff cot (3 



(7) 


takes the value Xt = \/Q (so called maximal mixing). For At = 1.5 TeV, Msusy = 1 TeV 
and mt = 175 GeV, one hnds Xt ~ 1.5, which yields lower Higgs masses than maximal 
mixing and so we are not typically choosing the most difficult scenarios that we could. 

Overall, for the squark and gaugino masses above, it might be that direct detection of 
the supersymmetric particles would not be possible and that the only new signal would be 
the detection of a Higgs boson. This makes the issue of whether or not at least one of the 
NMSSM Higgs bosons is guaranteed to be detectable at the LHG of vital importance. 

Finally, in our scan over parameter space, we restrict ourselves to the region mh± > 
155 GeV. For moderate tan/3, this means that t h^b and other possible charged Higgs 
signals would not reach the > 5a level at the LHG (see later discussion). 

^In addition, if the sparticles have large masses their contributions to the loop diagrams inducing Higgs 
boson production by gluon fusion and Higgs boson decay into 77 are negligible. 
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3 LHC prospects when Higgs-to-Higgs decays are for¬ 
bidden 

In this section, we update the no-lose theorem for NMSSM Higgs detection for the small 
portion of parameter space in which there are no decays of neutral Higgs bosons to other 
Higgs bosons. In the absence of such decays, the most relevant modes for detecting the 
neutral NMSSM Higgs bosons are those that have been earlier considered for the SM and 
for the MSSM [15-18]. These are (with £ = e,/x) 

2) associated Wh/a or tth/a production with in the final state; 

3) associated tih/a production with h/a ^ bb; 

4) associated bbh/a production with h/a ^ 

5) gg ^ h ^ ZZ^*^ 4 leptons; 

6) gg^h^ ITiyb) ^ i+i-uu; 

7) WW r+r"; 

8) WW WW^*\ 

In addition to these, we also include in our work the mode 

9) WW —h —>• invisible. 

LHC sensitivity to this mode has been studied by the ATLAS [16,19] and CMS [19-21] 
collaborations. In the present study, we employed the results of ref. [20] (fig. 25) that covers 
the Higgs mass range [100 GeV, 400 GeV] and includes systematic uncertainties. Details 
regarding how we treat all these modes are given in the Appendix. 

As regards the charged Higgs boson, it is well established (for early studies, see refs. [15, 
16]) that if f —>■ h^b decays are kinematically allowed then the will be relatively easily 
discovered in the decay of the top quark in tt events. These earlier studies have been updated 
by ATLAS in [22]. The resulting L = 300 fb“^ plot appears in fig. [TJ The least sensitivity 
to a charged Higgs boson occurs for tan ~ 4 10 where 5a is only attained for the MSSM 

parameter mA<135 GeV, corresponding to m/i±<155 GeV. Once m/i±>210 GeV, the Higgs- 
to-Higgs decay mode W^hi is typically kinematically allowed and in the NMSSM can 

have substantial branching ratio. Thus, in the NMSSM context the limits of fig. ^ on the 

in the region m/i±>210 GeV do not apply. 

We will study the complementarity between charged Higgs detection and neutral Higgs 
detection in the standard modes 1) - 9). We find that the smaller the lower limit on mh± for 
which we assume good significance for t ^ h^b detection, the smaller can be the minimum 
statistical significance for the neutral Higgs detection modes. 

In Ref. [4], a partial no-lose theorem for NMSSM Higgs boson discovery at the LHG 
(when Higgs-to-Higgs decays are forbidden) was established based on modes 1) - 8) above. 
There, we estimated the statistical significances {Nsd = S/\/B) for modes 1) - 8). For these 
results, it was especially critical that the tth with h ^ bb and IVIF-fusion modes [3) and 7), 
respectively] were included. Also important was mode 4), bbh with h . In the case of 

tth with h 65, we used the experimental study done by V. Drollinger at our request (see 
the Appendix) that extends results for this mode to Higgs masses as large as 150 GeV. The 
conclusion of ref. [4] was that, for an integrated luminosity of L = 300 fb“^ at the LHG, all 
the surviving points yielded Nsd > 10 after combining all modes. This means that NMSSM 
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Figure 1: We display the contours for 5a charged Higgs detection from [22], 


Higgs boson discovery by just one detector with L = 300 fb“^ is essentially guaranteed for 
those portions of parameter space for which Higgs boson decays to other Higgs bosons are 
kinematically forbidden. 

For the present paper, we have repeated the scan described above with the latest available 
LEP constraints (as incorporated in NMHDECAY [9] - see references therein) and with the 
latest ATLAS and CMS results for the discovery channels 1) - 9) (see the Appendix) and most 
recent 5cr curve for the as given in Eg. ^ For each Higgs state, we calculated all branching 
ratios using NMHDECAY. We then estimated the expected statistical signihcances at the 
LHC in all Higgs boson detection modes 1) - 9) by rescaling results for the SM Higgs boson 
and/or the MSSM h,H and/or A. The rescaling factors for the CP-even hi are determined 
by Ri, ti, bi = Tj, Qi and 7 j, the ratios of the VVhi, tthi, bbhi (or T~^T~hi), ggR and 77 / 1 * 
couplings, respectively, to those of a SM Higgs boson of the same mass. Of course \Ri\ < 1, 
but ti, bi, gi and 7 * can be larger, smaller or even differ in sign with respect to the SM. The 
reduced couplings for the CP-odd Higgs bosons (denoted by primes) are as follows: R[ = 0 
at tree-level; t'j and 6 ) are the ratios of the ^75 couplings for tt and bb, respectively, relative 
to SM-like strength. The quantities g'i and 7 ' are the ratios of the e x e' (e and e' being the 
polarizations of the gluons or photons) ggai or 770 * coupling strength to the e ■ e' gghsM 
or 'yyhsM coupling strength for = rrihsM- ^ detailed discussion of the procedures for 
rescaling SM and MSSM simulation results for the statistical significances in channels 1) - 9) 
is given in the Appendix. We will now summarize the results of this new no-Higgs-to-Higgs 
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scan. 


Only a few parameter choices are snch that decays of a nentral Higgs boson to any 
other Higgs boson are all forbidden when m^± > 155 GeV is reqnired. After restricting 
to snch parameters, it is extraordinarily difhcnlt to locate points that do not yield large 
statistical signihcance for LHC discovery of at least one Higgs boson while at the same 
time LEP constraints are not violated. We obtained a sample of 2455 points that had LHC 
signihcance (for L = 300 fb~^) in channels 1) - 9) below lOa. Most of these points had 
4< tan/3<10 — high tan/3 enhances prodnction cross sections for some of the Higgs bosons 
and will typically lead to visible signals. All points had LHC statistical signihcance above 
5 ct, thns establishing a no-lose theorem for points chosen consistent with LEP constraints 
and absence of Higgs-to-Higgs decays. Statistics on the important channels for these 2455 
points are snmmarized in tabled Note the importance of the channels 3), 4) and 7) for 
these most difhcnlt cases. 


Channel with highest S/\fB 

1 

2 

3 

4 

5 

6 

7 

8 

9 

No. of points 

0 

0 

343 

132 

0 

1 

1979 

0 

0 


Table 1: Most important channel for detecting the 2455 no-Higgs-to-Higgs-decays points 
that were most difhcnlt for LHC detection. 

In hg. H we give a scatter plot of the largest statistical signihcance achieved for a single 
nentral Higgs boson, A^^^, as a fnnction of mh± for the above 2455 points. (The absence 
of points with m/i±>210 GeV is dne to the fact that above this scale Higgs-to-Higgs decays, 
typically 02 —>■ Zhi and —>• hP^hi, wonld be allowed.) We see that the larger the valne 
of mh± for which a 5a signal can be established nsing tt prodnction with t ^ h^b decay, the 
larger the minimnm possible valne of the nentral Higgs bosons’ This means that the 

ATLAS and CMS gronps shonld work to maximize sensitivity to charged Higgs prodnction 
as well as to nentral Higgs prodnction. 

The point yielding the very lowest LHC statistical signihcance had the following param¬ 
eters, 

A = 0.0163; K =—0.0034; tan/3 = 5.7; 

/ieff = -284 GeV; Aa = -70 GeV; = -54 GeV, ( 8 ) 

which yielded m/j± ~ 155 GeV and nentral Higgs boson properties as given in table |21 Other 
points among the 2455 are similar in that the Higgs masses are closely spaced and below or at 
least not far above the WWjZZ decay thresholds, the CP-even Higgs bosons tend to share 
the WW/ZZ conpling strength (indicated by Ri in the table), conplings to bb of all Higgs 
bosons (the bi or 6 ' in the table) are not very enhanced, and conplings to gg and 77 (the gi 
or g[ and 7 ^ or 7 ' in the table) are snppressed relative to the SM Higgs strength. The most 
visible processes for this point had Nsd valnes above 6 . These were the WW h 2 ^ r+r“, 
WW —> ha —»• T^T~ and tth 2 ttbb channels. Overall, we have a qnite robnst LHC no-lose 
theorem for NMSSM parameters snch that LEP constraints are passed and Higgs-to-Higgs 
decays are not allowed. 

Another example point of possible interest is that giving the weakest signals when the 
charged Higgs mass is near the npper end of the spectrnm for which the does not decay 
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Figure 2: for the neutral Higgs bosons vs. mh± for m/j± > 155 GeV. Here, is 

the largest of the significances for any one of the neutral Higgs bosons in channels 1) - 9). 
We assume LHC luminosity of L = 300 fb“^. 


to other Higgs bosons. (This requirement is what restricts the upper range of mh± in hg. |21) 
The parameters for this point are: 

A = 0.0379; n = -0.0238; tan p = 8.25; 

/ieff = -119 GeV; Ha = -116 GeV; = -102 GeV , (9) 

yielding a charged Higgs mass of 216 GeV and neutral Higgs properties and statistical 
significances as listed in table El Note the strong signals in the WW hi _2 —>■ r+r“, 
tthi ^2 —ttbb and bbhs, 02 —> bbT~^T~ channels. 

These particular points in parameter space illustrate the general conclusion that it will be 
important that all the neutral Higgs detection modes that have been simulated by ATLAS 
and CMS really achieve their full L = 300 fb“^ potential. If the effective luminosity accu¬ 
mulated in modes 3), 4) and 7) were to all fall below L ~ 100 fb~ , then all single channel 
statistical significances for the most marginal points (as exemplified by the two tabulated 
points) would fall below 5 ct. Channel combination would be required to reach the > 5cr level. 
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Higgs 

hi 

h2 

hs 

Oi 

02 

Mass (GeV) 

99 

114 

145 

98 

134 

R 

0.49 

0.72 

-0.48 

- 

- 

ti or 

0.46 

0.65 

-0.64 

- 0.01 

0.18 

bi or 

1.71 

3.23 

4.49 

0.36 

5.59 

9i or 5 ' 

0.41 

0.56 

0.79 

0.02 

0.14 

7 * or 7 ' 

0.51 

0.75 

0.43 

0.01 

0.10 

BR{hi or ai —> hb) 

0.91 

0.90 

0.88 

0.92 

0.91 

BR{hi or Oi T^T~) 

0.08 

0.08 

0.09 

0.08 

0.09 

Chan. 1) S/^B 

0.00 

0.22 

0.20 

0.00 

0.00 

Chan. 2) S/y/B 

0.42 

0.80 

0.15 

0.42 

0.00 

Chan. 3) S/^/B 

3.52 

6.25 

5.39 

3.52 

5.39 

Chan. 4) S/VB 

0.73 

1.26 

3.86 

1.26 

3.86 

Chan. 5) S/^fB 

0.00 

0.15 

1.00 

- 

- 

Chan. 6 ) S/^ 

0.00 

0.00 

0.80 

- 

- 

Chan. 7) S/yfB 

0.00 

6.70 

6.54 

- 

- 

Chan. 8 ) S/^ 

0.00 

0.20 

0.25 

- 

- 

All-channel Sj^jB 

3.61 

9.29 

9.41 

3.76 

6.63 


Table 2: Properties of the neutral NMSSM Higgs bosons for the most difficult no- 
Higgs-to-Higgs-decays LHC point. In the table, Ri = ghiVV/QhsMVV, U = ghitt/9hsMth 
bi = g^bb!9hsMbb^ 9i = ghigg/ghsMjg and^y* = ghni/ghsMii ^Hsm = Similarly, t' and 
6 ' are the *75 couplings of to tt and bb normalized relative to the scalar tt and bb SM Higgs 
couplings and g[ and 7 ' are the a^gg and 0*77 e x e' couplings relative to the e • e' coupling 
of the SM Higgs. 


Higgs 

hi 

h2 

hs 

Oi 

02 

Mass (GeV) 

113. 

126. 

203. 

150. 

202 . 

R, 

0.75 

0.66 

0.01 

- 

- 

ti or t'i 

0.74 

0.65 

0.22 

- 0.01 

0.12 

bi or 6 ' 

1.38 

1.18 

- 8.11 

-0.78 

8.22 

9i or g'i 

0.72 

0.63 

0.37 

0.01 

0.08 

7 * or 7 ' 

0.76 

0.65 

0.93 

0.06 

0.05 

BR(hi or Qi bb) 

0.87 

0.80 

0.86 

0.91 

0.90 

BR(hi or Ui T+r“) 

0.08 

0.08 

0.09 

0.09 

0.10 

Chan. 1) S/VB 

1.91 

2.05 

0.00 

0.00 

0.00 

Chan. 2) S/VB 

4.99 

4.02 

0.00 

0.00 

0.00 

Chan. 3) S/VB 

7.84 

5.95 

0.00 

0.00 

0.00 

Chan. 4) S/VB 

0.20 

0.20 

7.40 

0.05 

7.40 

Chan. 5) S/VB 

1.33 

2.84 

0.31 

- 

- 

Chan. 6 ) S/VB 

0.00 

2.23 

0.19 

- 

- 

Chan. 7) S/VB 

6.67 

7.62 

0.00 

- 

- 

Chan. 8 ) S/VB 

1.14 

3.85 

0.00 

- 

- 

All-channel SjsjB 

11.73 

11.91 

7.41 

0.05 

7.40 


Table 3: Properties of the neutral NMSSM Higgs bosons for the second interesting point 
described in the text. Notation as in table |21 
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4 LHC prospects when Higgs-to-Higgs decays are al¬ 
lowed 

In this section we will consider the part of the parameter space complementary to that 
scanned in section IHl To be precise, we require that at least one of the following decay 
modes be kinematically allowed for some h and or a: 

i) h ^ h'h' , ii) h ^ aa , Hi) h , iv) h ^ aZ , 

v) h ^ , vi) a' ^ ha , vii) a hZ , viii) a . (10) 

(Recall that we do not consider —>■ W^h or —>■ W^a in defining the Higgs-to-Higgs de¬ 
cay parameter region.) The branching ratios for all these decays are computed by NMHDE- 
CAY. As in the previous section, we also allow for (but do not require) Higgs decays to 
gauginos. The large gaugino masses we employ imply that such decays are never important 
for the scans discussed in this paper. 

For most of these points it turns out that 5cr discovery of a neutral Higgs boson in at 
least one of the modes 1) - 9) is still possible. The number of parameter space points for 
which one or more of the decays i) — viii) is allowed, but 5cr discovery of a neutral Higgs 
boson in modes 1) - 9) is not possible, represents less than 1% of the physically acceptable 
points; in our scan we have found 3480 such points. In one sense, this small percentage is 
encouraging in that it implies that the standard LHC detection modes 1) - 9) suffice for 
most of randomly chosen parameter points. However, it should be noted that the fraction of 
points for which modes 1) - 9) suffice will decrease rapidly as the assumed LHC integrated 
luminosity is reduced. Further, the difficult parameter points are preferred on the basis 
of keeping fine-tuning modest in size [3]. (Modest fine-tuning means that is not very 
sensitive to GUT scale choices for the soft-SuSY-breaking parameters.) 

The parameters associated with these points for which all NMSSM Higgs bosons escape 
LEP detection and LHC detection in modes 1) - 9) occur throughout the broad range 
defined in eq. The scenarios associated with these points have some generic properties 
of considerable interest that make them worthy of further study. First, for all these 3480 
points, the and 02 are so heavy that they will only be detectable if a super high energy LC 
is eventually built so that Z ^ h^a 2 is possible, implying that LHC Higgs detection 

must rely on the lighter hi, ^2 and Oi states. The NMSSM parameter choices for which the 
latter cannot be detected at the LHC in the standard modes are such that there is a light, 
fairly SM-like CP-even Higgs boson (hi or h 2 ) that decays mainly to two lighter CP-odd 
or CP-even Higgs bosons (hi 2 —>■ OiOi or h 2 —>■ hihi). We will denote the parent SM-like 
CP-even Higgs boson by h^ and the daughter Higgs boson that appears in the pair decay 
by hL. 

We should discuss how it is that a light h^ will have escaped LEP detection. Consider 
the case of h^ = oi. First, sum rules require that the Zhiai {Zh 2 ai) coupling is small when 
the hi (h 2 ) coupling is near SM strength, implying that discovery in the e+e“ —>• Z* —>■ hiOi 
(e+e“ —> Z* —> h 2 ai) mode will not be possible. Second, e+e“ —> Z* —h 2 ai (hiOi) LEP 
constraints can be evaded in the NMSSM since the Oi can have sufficient singlet component 
that the Zh 2 ai (ZhiOi) coupling is small when the hi (h 2 ) is SM-like. For scenarios in which 
the h 2 is SM-like and decays primarily via h 2 —>■ hi hi, the hi is not observed at LEP because 
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of its weak ZZ coupling, while the ^2 mass is beyond the reach of LEP. II 

As we review the properties of specific bench mark points, it is useful to keep in mind the 
fact that detection of Higgs-pair final states at the LHC might be possible in certain cases. 
In particular, in [6, 7] we proposed that the LHC may be able to detect Higgs-pair hnal 
states using the WW —> > hihi —>• jjT^T~ production/decay mode. By this we mean 

that one of the final hi Higgs bosons is required to decay to a t^t~ pair where we identify 
the r’s through their leptonic decays to electrons and muons. The other final hi is required 
to decay to either bb (which we identify as jets) or, if bb is not a kinematically allowed decay, 
to jj or (with much smaller identification efficiency) where the r’s of this second pair 

would be tagged via their decays to single jets plus neutrinos. However, for a small fraction 
of the 3480 points, hn —^ cqcq decays are prominent but < 2mT-. For another small 
fraction, the hi has suppressed couplings to bb and r+r~. In either case, r triggering does 
not work and NMSSM Higgs detection at the LHC would probably be impossible. We will 
discuss this more in section El 

The distribution of the mass of the heavier SM-like Higgs {hn, where hn = hi or h 2 ) 
as compared to the mass of the lighter Higgs {hi, where hi = ai or hi) appearing in the 
hn hihi decay for the 3480 points from the scan described above is given in the top plot 
of fig. El We see that the SM-like parent Higgs mass uih^j lies in the [75 GeV, 155 GeV] 
interval while daughter Higgs masses range from near 0 up to close to mhfj/2. The 

middle plot of fig. Elshows that the parent hn always has fairly SM-like coupling iRh^l > 0-7 
to vector bosons; for most of the points, IRhhI is quite close to 1. The bottom plot of fig. El 
shows that for these points BR{hH —> hihi) is always substantial. The importance of the 
WW hn ^ hihi discovery mode is thus evident. 

Out of the above 3480 points, we have selected eight benchmark points, the properties 
of which are displayed in tables El and El that illustrate the cases where LHC detection of 
the NMSSM Higgs bosons in the standard modes 1) - 9) would not be possible. The first 
hve are such that the WW h^ ^ hihi —>• jjT^T~ detection mode might be effective. 
Points 6, 7 and 8 are chosen to illustrate cases where the hi appearing in the hnal state does 
not decay to either bb or r+r“, implying that the WW hn ^ hihi jjT~^T~ potential 
detection mode would not be useful. 

We now discuss in more detail the characteristics of these eight benchmark points. 

• Points 1, 2 and 3 are designed to illustrate hi —>• oiOi decay cases for a selection of 
possible hi and Oi masses. 

Point 1 is in the low-mass tail of the h^ mass distribution (see hg. El) at 90 GeV (al¬ 
though mhu as lows a 80 GeV is possible). For point 1, rria^ is below the bb threshold 
so that the main oi decay is to or jj. 

Point 2 and point 3 are at the two extremes of the central bulk of the h^ mass distri¬ 
bution of hg. El with nihfj = 100 GeV and 120 GeV, respectively. For these latter two 
points ruai is 20 GeV or 30 GeV; ai —>• bb and Ui —decays will be dominant 
and in the usual ratio. 

• Point 4 is such that the hi and /i 2 (with masses ruh^ = 97 GeV and ru/jj = 150 GeV) 

IIA similar situation arises in the case of a CP-violating MSSM Higgs sector [23]. There, the three Higgs 
bosons are mixed and parameter choices for which /12 ^ hihi decays are dominant can be found for which 
LEP constraints would not apply despite the fact that the hi is quite light. 
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hn^hjliL decays: mass pairings 
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share the WWjZZ coupling strength squared and both decay to aiOi. The ai decays 
to bb and in the usual ratio. Note that this point is an example for which ma^ is 
fairly large {nia^ = 45 GeV). 

• Point 5 illustrates a case in which it is the h 2 that is SM-like and it decays to hi hi. 
The hi —66 and hi —r+r“ decays are the dominant ones and are in the usual ratio. 
Although ruh-^ is rather small in this case, it would not have been seen at LEP due 
to its singlet nature. Nonetheless, BR{h 2 —^ hihi) is large due to the new trilinear 
NMSSM couplings. 

• For point 6, the hi is SM-like and decays via hi —> oiOi, but Oi —> 77 is dominant dne 
to the singlet nature of Oi. The 47 hnal state would provide a highly distinctive signal 
that should be easily seen at the LHC [14]. 

• Point 7 illustrates a case in which the h 2 is SM-like and decays via h 2 —hi hi. The 
new feature compared to point 5 is that the hi has rednced coupling to 66 and 

dne to the fact that parameters are snch that hi is almost entirely in nature [ 8 ]. 
** Obviously, the WW h 2 ^ mode would not be relevant for this type of 

scenario. We do not think that the resulting h 2 —4j signal could be isolated from 
backgrounds. 

• Point 8 illustrates a case in which the hi is SM-like and decays via hi —aiOi. It 
differs from earlier such points in that the ai is extremely light and decays mainly 
fo jj (j = s,c,g). Like for point 7, the WW —hi — jjT~^T~ detection channel 
would not be relevant. We would need to isolate a hi ^ Aj signal within a large QCD 
backgronnd. We do not believe this will be possible, especially given that each of the 
pairs of jets will have small mass and large boost, making separation of the two jets 
within each pair very problematical. 

This Oi would not have been seen at LEP in the hiOi mode for several reasons (for 
details see the references and discussions in [9]): 

1. First, the Zhiai coupling is very small because of the very SM-like nature of the 
hi. 

2 . Second, ma^ is below the threshold of any existing stndy of the ha type of mode 
at LEP. 

3. Third, the existing Zhi Zaiai limits (from OPAL) for the case where ai jj 
do not cover the mass regions corresponding to the valnes of either or maj. 
In fact, they extend only up to ~ 90 GeV and in any case do not apply when 
fUai is below the 2 GeV threshold for their searches. 

Finally, we note that axion searches do not apply since the Oi would not have been 
invisible in the detector — it decays promptly to visible jets. 


**A continuum of points of this type was discussed in ref. [9]. 
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Point Number 

1 

2 

3 

4 

5 

Bare Parameters 


A 

0.22 

0.4 

0.22 

0.67 

0.56 

K 

- 0.1 

-0.35 

0.59 

0.2 

0.1 

tan/3 

5. 

15. 

7.8 

4.1 

2.5 

/ieff (GeV) 

-520. 

-160. 

530. 

- 200 . 

-180. 

Aa (GeV) 

-580. 

-580. 

-920. 

-600. 

-440. 

A, (GeV) 

- 2.8 

-8.7 

- 2.1 

-30. 

172. 

CP-even Higgs Boson Masses and Couplings 


nih, (GeV) 

90. 

100 . 

119. 

97. 

40. 

Ri 

0.99 

0.97 

- 1.00 

0.69 

0.00 


0.99 

0.97 

- 1.00 

0.72 

0.05 

bi 

1.00 

0.90 

- 1.01 

0.31 

-0.35 

9i 

0.99 

0.97 

1.00 

0.74 

0.15 

7i 

0.99 

0.99 

1.00 

0.78 

0.11 

BRihi bb) 

0.08 

0.02 

0.01 

0.01 

0.93 

BRihi T+T“) 

0.01 

0.00 

0.00 

0.00 

0.07 

BR{hi aiOi) 

0.91 

0.97 

0.99 

0.99 

0.00 

To /12 (GeV) 

479. 

288. 

1431. 

150. 

125. 

i ?2 

0.16 

0.26 

0.00 

0.72 

- 1.00 

^2 

0.16 

0.26 

0.13 

0.70 

- 1.00 

^2 

0.19 

0.57 

-7.8 

1.10 

-1.03 

92 

0.16 

0.26 

0.12 

0.69 

0.99 

72 

0.19 

0.24 

0.08 

0.65 

0.99 

BR{h2 uioi) 

0.04 

0.44 

0.00 

0.97 

0.00 

BR{h2 hihi) 

0.25 

0.21 

0.00 

0.00 

0.92 

mhs (GeV) 

952. 

1016. 

2842. 

753. 

495. 

CP-odd Higgs Boson Masses and Couplings 


TOoi (GeV) 

10 . 

20 . 

31. 

45. 

144. 

A 

- 0.01 

0.00 

- 0.01 

- 0.02 

-0.06 

K 

- 0.22 

-0.85 

-0.53 

-0.40 

-0.40 

9'i 

0.15 

0.48 

0.19 

0.08 

0.06 

7( 

0.12 

0.11 

0.15 

0.49 

0.61 

BR{ai bb) 

0.00 

0.94 

0.93 

0.93 

0.85 

BRiai —!■ T+T“) 

0.83 

0.06 

0.07 

0.07 

0.08 

BR\ai jj) 

0.17 

0.00 

0.00 

0.00 

0.01 

TOaa (GeV) 

952. 

1018. 

1434. 

750. 

495. 

Charged Higgs Boson Mass 


TO/j± (GeV) 

954. 

1017. 

1432. 

742. 

487. 

LSP Mass 


m^o 

Xi 

453. 

136. 

476. 

113. 

82. 

Most Visible of the LHC Processes l)-9) 

2 (hi) 

5(/l2) 

2 (/ii) 

5(/l2) 

2 (/i 2 ) 

Nsd = S/VB of this process at L =300 fb“^ 

1.6 

0.7 

0.3 

0.5 

2.0 


Table 4: Properties of five scenarios for which LHC Higgs detection would only be possible in the WW 
hi_2 aiai jjT^T~ or WW h2 ^ hihi jjT~^T~ mode. The quantities Ri, U, bi, gi, 7 ^, t[, b[, g[ 
and 7 ' were defined in the caption of table |21 Important absolute branching ratios are displayed. Only the 
masses of the heavy / 13 , 02 and are given. The mass of the lightest neutralino (LSP) is also given. The 
second-to-last row gives the channel and Higgs boson yielding the largest Nsd = S/'/B in channels 1 ) - 9). 
The following row gives the corresponding Nsd for L = 300 fb“^. 
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Point Number 

6 

7 

8 

Bare Parameters 


A 

0.39 

0.5 

0.27 

K 

0.18 

-0.15 

0.15 

tanP 

3.5 

3.5 

2.9 

Meff 

-245. 

200 . 

-753. 


-230. 

780. 

312. 

A. 

-5. 

230. 

8.4 

CP-even Higgs Boson Masses and Couplings 


iTihi (GeV) 

94. 

57. 

95. 

Ri 

0.94 

-0.28 

1.00 

tl 

0.95 

-0.30 

0.99 

hi 

0.89 

0.01 

1.05 

91 

0.95 

0.33 

0.99 

7i 

0.96 

0.37 

1.00 

BR{hi jj) 

0.01 

0.93 

0.00 

BR(hi — > aitti) 

0.94 

0.00 

1.00 

ruh^ (GeV) 

239. 

125. 

483. 

R 2 

0.33 

-0.96 

-0.01 

t2 

0.30 

-0.95 

-0.36 

^2 

0.67 

-1.07 

2.84 

92 

0.29 

0.95 

0.37 

12 

0.30 

0.94 

0.68 

BR{h2 hihi) 

0.32 

0.93 

0.01 

mh 3 (GeV) 

562. 

731. 

821. 

CP-odd Higgs Boson Masses and Couplings 


mai (GeV) 

40. 

188. 

1 . 

t'l 

0.00 

0.04 

0.08 

b'l 

0.00 

0.53 

0.62 

9'i 

0.00 

0.04 

0.36 

li 

0.47 

0.31 

0.39 

BR{ai jj) 

0.00 

0.00 

0.95 

BR(ai — > /x/i) 

0.00 

0.00 

0.05 

BR{ai 77 ) 

0.98 

0.00 

0.00 

BR{ai XiXi) 

0.00 

0.99 

0.00 

TOaa (GeV) 

558. 

736. 

493. 

Charged Higgs Boson Mass 


mh± (GeV) 

560. 

727. 

485. 

LSP Mass 


TTlvO 

Xi 

211 . 

81. 

500. 

Most Visible of the LHC Processes l)-9) 

5(/l2) 

2 (/i2) 

5(/X3) 

Nsd = S/i/B of this process at L =300 fb“^ 

1.5 

1.3 

0.1 


Table 5: Properties of three representative scenarios for which LHC Higgs detection would not even be 
possible in the WW ^ oiQi ^ jjT^T~ or WW ^ /12 ^ ^1^1 ^ jjT~^T~ modes. 
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5 Collider Implications 

In the previous section, we have established the probable importance and possible necessity 
of detecting a fairly SM-like, relatively light Higgs boson, hu, in a Higgs-pair decay mode, 
hn —^ hihi- In this section, we discuss possible ways in which such detection might be 
possible at various different kinds of colliders, with some emphasis on the LHC. The best 
means for such detection at hadron colliders will depend strongly upon the hi decay channels. 
Detection of the hn in the Higgs-pair hnal state at an e+e“ or 77 collider will be less 
dependent upon precisely how the hi decays. 

The LHC 

At the LHC it will presumably be highly advantageous to use WW fusion production 
for the hn- Not only is the associated production cross section quite competitive with other 
production mechanisms for hn in the G [80 GeV, 150 GeV] mass range of relevance, 
but also the ability to tag the spectator jets will certainly make backgrounds much more 
manageable. 

When > 2mb, we advocate employing the h^hL ^ bbT~^T~ hnal state. In this hnal 
state an approximate mass for the hu can be computed using the visible particles in the hnal 
state to compute an ehective mass, M^jv+r- • In addition, restrictions on the visible mass of 
each hi can also be imposed. (In the analysis, one will need to choose a hypothetical value 
for and examine the Mjj^+T— mass distribution for a peak. This process will have to be 
repeated for all possible values. Only the choice with near the actual value would 
reveal a peak in For the hnal state, the two main backgrounds appear to 

be: (i) gg ^ tt ^ bbW~^W~ bbT~^T~ -|- ^t, in association with forward and backward jet 
radiation and (ii) Drell-Yan -|- jets production. Monte Carlo simulations performed 

to date show that 6 -tagging does not seem to be necessary to overcome the a priori large 
Drell-Yan + jets background. It is eliminated by stringent cuts for hnding the highly 
energetic forward / backward jets characteristic of the WW fusion process. To the extent 
that the main background will then come from tt production, it is not useful to specihcally 
6 -tag the 6 jets since the tt background will also contain 6 jets. Thus, it is appropriate to 
focus on a generic WW jjT^T~ -|- hnal state. Such a hnal state can be experimentally 
isolated with high efficiency by identifying two r’s from one Oi using the leptonic decay 
modes for both r’s while requiring two ( 6 ) jets from the other Oi. All these particles should 
be required to be quite central. 

If ruh^ < 2mb, but above 2mr, the dominant hnal state is . However, 

an ehective hn mass is very difficult to reconstruct in this channel. Previous work suggests 
that it will be best to employ the hihi ^ hnal state (which typically has a small but 

usable branching ratio) where one of the h^’s decays to jj = cc, ss or gg. This is extracted 
experimentally by again identifying two r’s in their leptonic decay modes and two jets. 

Preliminary simulations for the jjr^r“ signal have appeared in [ 6 ,7] for a few represen¬ 
tative benchmark points (diherent from those appearing in the tables of this paper). Aside 
from imposing stringent forward / backward jet tagging cuts to eliminate the Drell-Yan 
r+r“ + jets background, it was required that the two additional jets (from one of the h^s) 
and the two opposite sign central leptons {i = e, g) coming from the the r’''r“ emerging from 
the decay of the other all be quite central. Additional observations are the following: 

• In the case of > 2mb, the jj mass is fairly high, efficiencies for identifying such j’s 
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were found to be high and the momenta of each j were relatively well determined. The 
latter implies that the r+r” mass can be reconstructed by assuming that all missing 
transverse momentum is to be associated with the neutrinos in the r+r“ —> + u's 

decays. Efficiencies for the overall reconstruction of this kind of event are therefore 
reasonably high. 

• In the case of 2771,- < < 2mh, the primary source of the jj is from cc + ss + gg 

decays of one of the h^s (neglecting the inefficiently tagged and poorly reconstructed 
contribution coming from with two r —> j' + IpT decays, where the j is 

a single pion or similar hadronic resonance). However, the jj pair mass is quite low 
(< 2mh) and separate identihcation of the two jets was found to be rather inefficient. 

For this case, the relevant hnal state branching ratio is BR{hihi jjT^T~) = 2 x 
BRiJiL —> cc + ss + gg) BRihi —> which is similar in size (typically) for 

2 m,- < mh^ < 2mb to what one hnds when uih^ > 2mb for BRiJiihL —^ bbT~^T~). 

We reiterate that since the hi will not have been detected previously, we must assume a 
value for m^,^ to perform the analysis. We then look among the central jets for the combi¬ 
nation with invariant mass Mjj closest to m^,^ (no 6-tagging is enforced, 6’s are identihed as 
non-forward/backward jets). We then compute the MjjT-+T-~ invariant mass using the four 
reconstructed four-momenta for the two j’s and two r’s and look for a bump in the distribu¬ 
tion. This process is repeated for densely spaced m/,^ values and we look for the rrih^ choice 
that produces the best signal. 

In our earlier simulations of points with m/,^ G [115 GeV, 130 GeV], the typical result 
found (for the assumed ruh^ chosen to agree with the actual was a sizable bump in 

MjjT+r- coming from the signal in the MjjT+r- range from roughly 40 GeV to 130 GeV. The 
ti background produces a huge peak in Mjjr+r- at high mass with a rapidly falling tail in the 
MjjT+r- ^120 GeV region. The crucial issue is the precise shape of this tail and exactly how it 
extends into the low-Mjj^+T— region where the signal peak resides. The early results of [6,7], 
based on UAl detector resolutions and efficiencies, found that the low-MjjT+r- tail from tt 
does not overlap signihcantly the signal bump, which (for L = 300 fb~^) typically contained 
between 500 and 2000 events in Mjj^+T— G [40 GeV, 130 GeVj. More recently, members of 
the ATLAS collaboration [24] have examined the MjjT+r- signal using the ATLAS simulation 
programs and somewhat different cuts. They hnd that the ti background extends over the full 
range where the signal resides. This question is now being examined using full simulations 
in the context of the ATLAS detector and improved cuts [24,25]. We are unable at this 
time to say whether or not the signal will emerge above the background in a statistically 
signihcant and reliable way. 

In principle, one could explore hnal states other than the hihi —> jjT^T~ mode. How¬ 
ever, all other channels will be much more problematical at the LHG. A 46-signal would be 
present for ruhi^ > 2mb but would be burdened by a large QGD background even after imple¬ 
menting 6-tagging. Meanwhile, for 2mr < < 2mb the 4r-channel would typically have 

a large branching ratio and we could look for it in the mode where all r’s decay leptonically. 
However, in this mode it would not be possible to reconstruct the hn resonance mass and 
backgrounds would be large. 

It should be clear that without the ability to tag one r+r“ pair as part of the hn 
reconstruction process, the background would be much larger. Thus, it is only for mh^ > 2mr 
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that there is a possibility to isolate the WW ^ hn ^ hihi signal. We are very pessimistic 
regarding isolating a signihcant signal in a 4j hnal state as appropriate when < 2mr or 
in those special cases where > 2mr but hi —^ jj (j = c, s, g) decays are dominant. 

The Tevatron 

At the Tevatron, WW —> hn production has a rather low cross section. Only the gg hu 
cross section is sizable (of order 1 pb for a SM-like hu with ~ 100 GeV). For the case 
of > 2mfe, one would again employ the hihi —^ bbT~^T~ final state. However, forward / 
backward jet tagging could no longer be used to reduce the Drell-Yan +jets background 
without also severely affecting the gg fusion signal. The best means for discriminating against 
the DY background would probably be to use 6-tagging. Of course, this will not reduce the 
dominant tt background relative to the Higgs signal. Detailed simulations will be required 
to see if a signal can be extracted.. A group [26] including CDF experimentalists is working 
on simulating typical cases with 2mr < nih^ < 2mb. They are currently focusing on the 
t+t-t+t~ hnal state where one r is identihed through its r —p -|- $t decay and all other 
r’s are identihed using either an isolated hadronic decay signature or a lepton decay. 

An e'*~e~ linear collider 

At an e^e“ collider, it will be possible to detect any relatively light Higgs boson with 
substantial ZZ coupling using the e’''e“ Z* ^ ZX hnal state and searching for the promi¬ 
nent peak in Mx that would arise if a Higgs boson is present. This technique is completely 
independent of the Higgs decay mode. Once a peak is found, it would be straightforward to 
isolate the hn —^ hihi hnal state in various hi decay modes and check if the hi branching 
ratios are consistent with expectations for a light hi or oi of the observed mass. In addition, 
^hLi QzzhH) ghnhLhL be measured with considerable precision. This would 

allow precision tests of the NMSSM model structure, especially if part of the supersymmetric 
particle spectrum is also accessible. 

A 77 collider 

Another facility of particular interest for the kind of scenario presented here will be a 77 
collider. Since the hn is typically quite SM-like, it will have a very substantial production 
rate in 77 collisions. A recent study [27] shows that a very substantial signal for the 77 — 
hn —^ hihi process will be present above a very small background (after appropriate simple 
cuts) in the main hihi ^ bbbb and hihi bbr'^T~ hnal states. Excellent determinations 
of both nihfj and will be possible and the 77 coupling of the hn will be very precisely 
determined, as will the hn —^ hihi coupling strength. 
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6 Conclusions 


In summary, we have explored the NMSSM model parameter space, looking for Higgs sector 
scenarios consistent with LEP exclusions that might be unexpectedly difficult to probe at the 
LHC in the conventional modes that have been explored for the SM and the MSSM. We have 
found that generic points in NMSSM parameter space are such that Higgs-to-Higgs decays 
are present. This is a crucial issue since hadron collider signals for Higgs bosons decaying 
to other Higgs bosons will typically be much more difficult to extract in the presence of 
backgrounds than signals in the conventional modes studied for the SM/MSSM scenarios. 

In section El we considered NMSSM parameter points for which decays of neutral Higgs 
bosons to other Higgs bosons were not present. For this small fraction of parameter space, 
we are able to show that the conventional SM/MSSM Higgs boson discovery modes 1) - 
9) (as listed at the beginning of the section Ej) are sufficient (assuming L = 300 fb“^) to 
guarantee that at least one NMSSM Higgs boson will be detectable at the > 5a level at the 
LHC. The worst point yielded signals between 6a and 7a in several of the standard modes for 
several different Higgs bosons. However, the limited statistical signihcance for these signals 
means that if the effective integrated luminosity falls below L ~ 100 fb~^ or if backgrounds 
are larger than found in the simulations, then this ‘no-lose’ theorem would fail. In any case, 
for the most difficult no-Higgs-to-Higgs points Higgs discovery will not be easy or quick — 
considerable thoroughness and patience will be required. The interplay between different 
detection channels and different Higgs states will be crucial. Good statistical significance 
might only be achieved by combining a number of channels. 

The vast bulk of physically acceptable NMSSM parameter choices are such that Higgs-to- 
Higgs decays are present. The focus of this paper has been to isolate those cases where these 
decays reduce statistical signihcances in all the standard modes to a level such that there is 
no > 5a level signal in any standard mode for any Higgs boson. In section EJ we presented 
eight sample points for which all the standard modes had very low statistical signihcance 
and detection of a SM-like hu decaying to a pair of lighter h^’s would provide the only 
possible signal. Five of the sample points were such that the hnal state of interest would be 
hihi —*■ jjT~^T~ (for 2mr < mh^ < 2mf,) or hihi —> bbT~^T~ (for 2mi, < rrih^)- We noted 
the potential importance of the LHC channel WW hu ^ for detecting 

a Higgs signal in these cases. Two of the other three points were such that only > jj 
(j = c,s,g) decays were present and the last point was such that —>• 77 . In the latter 
case, the 47 hnal state would provide a very clean LHC signal. In the former two cases, we 
are unable to envisage a technique for discovering any of the Higgs bosons at the LHC. 

In section El we pursued further the issue of Higgs detection at the LHC for cases like 
those of the hrst hve sample points noted above. We discussed the nature of and techniques 
for extracting the WW hn —>■ hihi ^ LHC signal. As noted there, this signal 

is being actively worked on in collaboration with members of ATLAS. We also presented a 
brief summary of the Higgs boson signal at the Tevatron based on gg —> hfj fusion. Finally, 
we summarized why it is that at an e+e“ collider or 77 collider it will be far easier to detect 
hfj production followed hj hn ^ decay than at a hadron collider. For instance, at the 
ILC, discovery of a light SM-like h is guaranteed to be possible in the Zhn hnal state using 
the decay-independent recoil mass technique [28]. 

Regarding the scenarios for which only the WW hn —^ hihi JjA+r" channel 
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might provide a signal at the LHC, we note that the main issue will be whether the back¬ 
ground from ti production (which we believe is the primary background after appropriate 
cuts requiring highly energetic forward / backward jets to eliminate the DY + jets 

background) will extend to low values of the reconstructed MjjT+r- mass where the signal 
resides. To answer this question requires a very full simulation. However, it is essential 
that the ATLAS and CMS groups attack this problem vigorously since, in the worst case 
scenarios, this signal will be the only evidence for Higgs bosons at the LHC. Once the LHC 
is operating, the tt background can be more completely modeled and the signihcance of any 
enhancement observed in the MjjT-+T— distribution more reliably assessed. However, even if 
a fully trustworthy signal is seen at the LHC, a future ILC will probably be essential in order 
to conhrm that the enhancement seen at the LHC really does correspond to a Higgs boson. 

We should also note that, for parameter space points of the type we have discussed here, 
detection of any of the other NMSSM Higgs bosons is likely to be impossible at the LHC and 
is likely to require an ILC with above the relevant thresholds for h'a! production, 

where h' and a' are heavy CP-even and CP-odd Higgs bosons, respectively. 

Although the scan results presented here were done for sparticles (except possibly the Xi) 
that are fairly heavy, we do not believe the results will change signihcantly if the sparticles 
are as low in mass as current LEP and Tevatron bounds. This is because the primary 
issue is how the SM-like Higgs boson (which must have mass below roughly 150 GeV when 
perturbativity up to the GUT scale is imposed) decays. Its decays will not be signihcantly 
affected by sparticles with masses even slightly above current limits. 

At the LHC, if SUSY is discovered and WW WW scattering is found to be pertur¬ 
bative at WW energies of 1 TeV (and higher), and yet no Higgs bosons are detected in 
the standard modes, a careful search for the signal we have considered should have a high 
priority. 

Finally, we should remark that the hn —^ hihi search channel considered here in the 
NMSSM framework is also highly relevant for a general two-Higgs-doublet model, 2HDM. It 
is really quite possible that the most SM-like CP-even Higgs boson of a 2HDM will decay 
primarily to two CP-odd states. This is possible even if the CP-even state is quite heavy, 
unlike the NMSSM cases considered here. If CP violation is introduced in the Higgs sector, 
either at tree-level or as a result of one-loop corrections, then hn —^ hihL decays will 
generally be present (as, for example, in the CP-violating MSSM [23]). The critical signal 
will be the same as that considered here. 
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Appendix A: Summary of ATLAS and CMS simulations 
employed and rescaling procedures 

We had a large number of experimental simulations available for each of the standard dis¬ 
covery channels 1) - 9). Because of the need to go to L = 300 fb~^ in order to achieve a 
hrm no-lose theorem for NMSSM Higgs discovery in the absence of Higgs-to-Higgs decays, 
whenever available we employed results for CMS or ATLAS for L = 100 fb~^ rather than 
low luminosity, L = 30 fb“^, results. In some channels, the CMS results indicated greater 
discovery potential than ATLAS results and vice versa. We always employed the best single 
detector results. (That is, we do not double the statistics assuming two detectors.) We did 
not make use of any studies other than those performed by the ATLAS and CMS detector 
collaborations. We do not attempt to give all the different simulations considered but only 
summarize those we actually used for each of the nine standard channels. We apologize in 
advance for not referencing all the experimental (and theoretical) studies that we did not 
end up using. 

To be conservative, we always employed results obtained for the case where the radiative 
correction “iL factors” for the signal and background were unity: Ks = 1. and iCs = 1. At 
the LHC, it is almost always the case that the actual K factors for the signal and background 
(before cuts) for a given channel are such that ^ improves upon their inclusion. But, using 
the K factors obtained before cuts is unreliable since the K factors can easily be sensitive 
to the cuts and selection procedures employed by the experimental groups. Eventually, full, 
process-specihc Monte Carlos will be available at NLO that will allow K factor evaluation 
after cuts, at which time this kind of study could be repeated in order to see if the radiative 
corrections have signihcant impact. For a recent summary of LHC radiative corrections 
related to Higgs production and decay, see [29] and references therein. 


Table 6: Resolutions for combining overlapping signals in a given channel. 


Channel Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

m 

0.01 

0.01 

0.10 

0.15 

0.01 

0.01 

0.10 

0.10 

0.10 


Finally, we must account for the fact that the different hi and a* can have a range of 
different masses, sometimes overlapping, sometimes not. Thus, signals in a given discovery 
channel from different scalars and/or pseudo-scalars can overlap within the experimental 
resolution. In this case, the overlapping signals should be combined. We have chosen to 
combine the scalar and/or pseudo-scalar signals at different masses following the procedure 
of ref. [30], section 5.4, using a channel-dependent resolution. In particular, we have chosen to 
employ (in the notation of [30]) CFm{i) = 2A(i) xm {m being the Higgs mass, and i denoting 
the channel) with the A(i) values as given in table [Ul A particularly relevant example is 
channel 4) (in the sense that there is often overlap between scalar and pseudo-scalar Higgs 
boson resonance signals which individually have a useful level of signihcance). For channel 
4) we estimated A = 0.15 from £g. 19-61 in [16] at high luminosity and extrapolated to 
niA < 150 GeV. 
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Channel 1); For ^ h —> 77 we employ L = 100 fb“^ resnlts analogons to those for 
L = 30 fb“^ contained in hg. 1 of [31]. For onr purpose it is crucial to avoid summing 
over the h —77 and Wh + tth —>■ channels. This is because production 

rates in these two channels are scaled differently in the NMSSM, the hrst being scaled 
by the factor gf and the second by a combination of Rj and tj. We thank R. Kinnunen, 
K. Lassila-Perini and A. Nikitenko for providing us with this separation in a series of 
email communications. For L = 100 fb~^, the resulting S/\/B values for the gg —^ hsM 
fusion process alone are summarized in table[Z|below for the assumption that K factors 
for signal and background are both unity: Ks = Kb = I- 


The production rates in the gg channel must be corrected for non-SM-like gg couplings 
of the Higgs bosons. We must also account for differences in the 77 branching ratio 
relative to that of the SM Higgs boson. Thus, the tabulated entries are to be multiplied 
by gfBR{hi 'y'j)/BR{hsM 77 ) for the hi^ 2,3 and by g['^BR{ai 'j'y)/BR{hsM 
77 ) for the ai, 2 - The L = 300 fb~^ results are obtained by scaling the results so obtained 

, Ton. 1/2 

by 


300 fb" 

100 fb 




Table 7: gg ^ h ^ jr CMS, L = 100 fb"\ Ks = Kb = I 


m [GeV] 

100 

110 

120 

130 

140 

150 

S/^ 

4.2 

6.0 

6.8 

8.2 

7.0 

5.2 


Channel 2): For Wh + tth —> 77 , we employ the CMS L = 100 fb“^ results for this 
separate channel, as provided to us by R. Kinnunen, K. Lassila-Perini and A. Nikitenko. 
These are tabulated in table |H1 Since the S/\fB in these SM-Higgs simulations came 
about 50% from the Wh channel and about 50% from the tth channel, we rescale the 
production rate for this process by | {Rf -|- tf) for hi^ 2,3 or for the 01 , 2 - Including 
the correction for the 77 branching ratio, the tabulated results are rescaled by | -|- 

t^)BR{hi 'y'y)/BR{hsM ^ 77 ) for the hi^ 2,3 and by \ t[^BR{ai 'y'j)/BR{hsM ^ 
77 ) for the ai^ 2 - The L = 300 fb“^ results are obtained by scaling the results so 


obtained by 


300 fb- 


100 fb 




Table 8 : Wh/tih -f-fl: CMS, L = 100 fb“\ Kg = I, Kb = I 


m [GeV] 

80 

90 

100 

110 

120 

130 

140 

150 

s/Vb 

9.4 

10.6 

10.9 

14.8 

15.7 

13.2 

10.4 

8.2 


Channel 3): For tth —> ttbb, we employed results supplied by V. Drollinger based on 
extension of the work in ref. [32] to the much larger Higgs mass range required for our 
NMSSM study. We are very grateful for these additional results, which were absolutely 
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critical to our study, and for the collaboration of V. Drollinger in checking the hnal 
table El below, including: the extrapolation to L = 100 fb“^; the change from Kb = 1-9 
to Kb = ^ (the standard we employ in this paper); and the removal of the SM result 
for BR{hsM bb) — the results of table El are to be multiplied by BR{h, a —>■ bb) 
and not the ratio to the SM Higgs branching ratio. V. Drollinger emphasizes that 
the extrapolation to L = 100 fb“^ has ignored beam pile-up which might cause some 
diminution in 6-tagging efficiency at the higher L = 100 fb~^ luminosity. (This will be 
studied during preparation of the CMS TDR.) Thus, we have been somewhat cautious 
in extrapolating table El fo the full L = 300 fb~^ luminosity by employing the factor 

As regards rescaling this table for the various NMSSM Higgs bosons, the 

bb) 


300 fb- 


1/4 


100 fb^ 

resnlts given are to be mnltiplied by tfBR{hi bb) for hi, 2,3 and by t'^BR{ai 
for oi 2 - 


Table 9: tth —^ ttbb\ L = 100 fb , Kg = Kb = 1, quoted for BR{h bb) = 1 


m [GeV] 

80 

90 

100 

no 

120 

130 

140 

150 

siVb 

17.9 

15.0 

14.1 

12.3 

12.7 

13.7 

11.3 

10.6 


Channel 4): For bbh/a —>■ bbT~^T~ we have employed the experimental studies presented 
in [16] (as contained in the L = 100 fb“^ curve of £g. 19-62 and also using information 
in tables 19.35/36). These results were repeated (for the mass range below 500 GeV 
where we employ them) in the Les Honches workshop stndy of [19] £g. E.15. The 
estimation of the statistical signihcances using £g. 19-62 of [16] for this channel reqnires 
considerable discussion. 

Figure 19-62 gives the 5a contours in the tan/? - tua plane of the MSSM. The critical 
issue is what fraction of these 5a signals derives from gg ^ H/gg ^ A production and 
what fraction from associated bbH/bbA production, and how each of the gg fusion and 
bb associated prodnction processes are divided up between H and A. For the former, 
we turn to table 19.35 of ref. [16]. There, we see that it is for cuts designed to single out 
the associated production processes that large statistical signihcance can be achieved 
and that such cuts provide 90% of the net statistical signihcance of Ngn = 8.9 (3.9 for 
gg fusion cuts combined in quadrature with 8.0 for bbH + bbA associated production 
cuts) for rriA = 150 GeV and L = 30 fb“^. (For the associated production cuts, the 
table of [16] shows that the contribution of the gg fusion processes to the signal is very 
small.) The percentage of Ngn deriving from (/-fusion cuts is even smaller at high 
tua- For mn ~ rnA € [100 GeV, 500 GeV], a conservative choice is then that 90% of 
the statistical significance along the contours of £g. 19-62 comes from the associated 
production cut analysis. With this choice, the 5 ct contour at L = 100 fb“^ from fig. 
19-62 of ref. [16] corresponds to a 4.5 ct contour for associated bbH -|- bbA production 
alone. Since the values of tan /? along this contour are large, we can separate the H 
and A signals from one another by using the following properties of the MSSM within 
which £g. 19-62 of ref. [16] was generated: (a) BR{H r+r“) ~ BR{A — t~^t~) ~ 
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0.09; (b) the bbA and bbH couplings are very nearly equal and scale as tan/?; and (c) 
rriA ~ mn within the mass resolution. As a result, the net signal rate along this 
contour is approximately twice that for bbA or bbH alone. Thus, Nsd = 2.25 would 
be achieved for bbA or bbH along this contour were ttia and mn widely separated. 
Dehning the value of tan/? as a function of tua shown by the 100 fb“^ curve of £g. 
19-62 in ref. [16] as tan/?2.25(^7 .a), we compute 


A^5£)(tan/? = 1, m) 


2.25 


tan/?2,25(ni) 


( 11 ) 


These are the numbers tabulated in tableElbelow (where, for convenience, we include 
an extra factor of 100). 

The above procedure is conservative in that it assumes no contribution to the 
channel Nsd from the gg fusion processes. We have not attempted to include the latter 
production process, since the mode is only useful in finding 5 ct contours when 

the bb Higgs coupling is highly enhanced, in which case the gg fusion process will make 
a relatively very small contribution. 


Finally, the results of ref. [16] assumed Ks = Kb = ^ and assumed the MSSM mass- 
independent value for the branching ratio, BR{H,A —> = 0.09. Putting all 

this together, the results of table HU] must be rescaled by the factor O.Olb'fBR^hi —> 
r+r“)/0.09 for hi 2,3 and by 0.01b'j^BR{ai —>• r’''r“)/0.09 for 01^2- The L = 300 fb“^ 

ronnru-llV2 

results are obtained by scaling the results so obtained by 


300 fb- 

100 fb 




Table 10: bbh or bba bbrr at g^jh^^a ~ 9bbH = !)• T = 100 fb Ks = Kb = 1 


m [GeV] 

100 

110 

120 

130 

140 

150 

200 

S/y/B (xl02) 

3.7 

4.2 

4.4 

4.5 

4.7 

4.6 

3.1 

m [GeV] 

250 

300 

350 

400 

450 

500 


S/y/B (xl02) 

2.1 

1.3 

1.0 

0.8 

0.7 

0.6 



Channel 5: For gg ^ h ^ > Ai^Uuiy, we employ the CMS “no K-factors”, L = 

100 fb“^ plot supplied to us by R. Kinnunen. For Higgs mass below 500 GeV, only the 
M mode is present on the plot. For masses from 500 GeV up to 1 TeV, the tabulated 
numbers were obtained by combining in quadrature the plotted results for the Ai and 
2i2u modes. See also, [15]. (The GMS L = 30 fb“^ results appear in fig. 1 of [31] and 
figs. 12 and 13 of [20].) The results quoted in table [TT] assume Ks = Kb = 1. The 
tabulated values are to be multiplied by RfBR{hi —> ZZ^*'>)/BR{hsM —^ ZZ^*'^) for 
hi,2,3. We assume no contribution from this mode to 01,2 corresponding to the absence 

r -iGA 

of tree-level ZZoi couplings. For L = 300 fb“^, we scale by the factor 

Ghannel 6): For gg ^ h —>■ WW^*'^ —> we again employ the GMS Ks = Kb = I, 

L = 100 fb“^ plot supplied to us by R. Kinnunen. The Huu signal is the only one 
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Table 11: gg ^ h ^ ^ M, iiuu: L = 100 fb-\ Ks = Kb = I 


m [GeV] 

100 

120 

130 

140 

150 

160 

170 

180 

190 

200 

s/Vb 

2.7 

5.3 

13.2 

22.1 

27.8 

9.4 

5.5 

20.7 

25.1 

26.1 

m [GeV] 

250 

275 

350 

400 

500 

600 

700 

800 

1000 


S/^ 

21.6 

17.6 

22.7 

21.6 

21.5 

17.1 

13.6 

11.1 

9.3 



present for [120 GeV, 250 GeV]. At mass=300 GeV, both Huiy and £ijj are present, 
and we combine them in qnadratnre. For the masses of 600 GeV and 800 GeV, only the 
£ijj signal is present. The resnlts that we obtain in this way from the GMS plot areas 
tabulated in table IT^ below. (The GMS L = 30 fb“^ results appear in £g. 1 of [31] and 
hgs. 12 and 13 of [20].) For NMSSM Higgs statistical signihcances at L = 100 fb~^, 
we multiply the values in the table by gfBR{hi —>• WW^*^)/BR{hsM WW^*^) for 
the hi^ 2 , 3 - This channel is absent at tree-level for the ai, 2 - In going to L = 300 fb~^, 

r -iG /2 

results obtained in this way were multiplied by ^_i 


Table 12: gg ^ h ^ WW^*^ ££uu, £ujj: L = 100 fb-\ Ks = Kb = I 


m [GeV] 

120 

130 

140 

150 

160 

170 

180 

s/Vb 

5.1 

9.8 

17.8 

21.9 

47.0 

34.4 

24.1 

m [GeV] 

190 

200 

250 

300 

600 

800 


S/^fB 

19.5 

16.9 

7.9 

19.4 

14.2 

11.3 



Ghannel 7): For the W~^W 
ATLAS results of 


■> h ^ T~^r~ channel, we employed the table 10, L = 30 fb“ 

r -iM/2 

, rescaled to L = 100 fb“^ by the factor of 


100 fb~ 
30 fb- 


The 


results of this rescaling are given in tableUHlbelow. The tabulated values are multiplied 
by R^BR{hi —> t^t~)/BRQism t^t~) for the NMSSM hi^ 2 , 3 - This mode is not 
present (at tree-level) for the ai^ 2 - In going to L = 300 fb“\ we multiplied the results 

so obtained by the somewhat conservative factor of 


300 fb " 

100 fb 




Table 13: WW r+r”: L = 100 fb-\ Ks = Kb = 1 


m [GeV] 

110 

120 

130 

140 

150 

S/^ 

6.7 

10.4 

10.4 

8.7 

4.4 
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Channel 8): For WW h —>■ WW, we employed the results in table 7 of [33] in the last row 
labeled “combined statistical signihcance”. These results were those obtained for L = 
10 fb“^. Since the main final state contributors to the statistical significances given for 
L = 10 fb“^ were the WW —*• ee, /i/r and e/r hnal states, we felt that these results could 


safely be scaled up to L = 100 fb ^ using the factor 


100 fb- 


-ll V2 


10 fb 




The results of this 


scaling are tabulated in table^^ In addition, there was a specialized neural net analysis 
for the limited mass range of [115 GeV, 130 GeV] [34]. The results corresponding to 
L = 100 fb“^ from table 5 of this analysis are given in the parentheses in table El In 
the [115 GeV, 130 GeV] mass range, we have employed the (stronger) neural net result. 
Entries in table HU are to be multiplied by R}BR{hi —> W^W~) / BR{hsM —*■ W^W~) 
for the hi^ 2 , 3 - The process is absent at tree-level for the ai^ 2 - In going to L = 300 fb“^, 

we have been somewhat conservative and scaled the results so obtained by the factor 

1 1 1/4 
300 fb-i 


100 fb“ 


Table 14: WW WW: L = 100 fb“\ Ks = Kb = I 


m [GeV] 

110 

115 

120 

125 

130 

140 

s/Vb 

2.5 

(5.6) 

6.6 (9.7) 

15.7 

13.9 (20.5) 

18.6 

m [GeV] 

150 

160 

170 

180 

190 


S/y/B 

26.5 

34.8 

34.8 

27.8 

21.5 



Cha nn el 9): For invisibly decaying Higgs bosons, there are two experimental studies. The 
hrst is that of [21] covering the Higgs mass range [100 GeV, 140 GeV]. This study was 
recently extended to a larger Higgs mass range in [20]. Both studies were performed 
for L = 10 fb“^. Since we are uncertain that these results can be easily employed at 
higher L, our program currently assumes that only L = 10 fb“^ of data is accumulated 
for the WW h —> invisible mode. The appropriate procedure for the results quoted 
in [20] as based on [19] is as follows. The raw S/\/B of this latter reference agrees 
(at Higgs mass = 120 GeV) with that in [21]. However, [20] includes a systematic 
3% uncertainty in the background and computes S/ -y/HTiT((pQ 3 ^j 2 obtain the 95% 
GL limits of their fig. 25. This we believe is the more reliable way of estimating the 
significance of the signal given the amorphous nature of the backgrounds. In table El 
we give the significances after accounting for the background systematic uncertainty. 
These are extracted from fig. 25 of [20] using the formula S/\/B + (0.035)2 ^ 1.96/^^, 
where = BR{h —>• invisible )is the quantity plotted. 























Table 15: WW h ^ invisible: L = 10 fb \ Ks = Kb = 1, BR{h —> invisible) = 1, SM 
WWh coupling 


m [GeV] 

120 

150 

200 

250 

300 

350 

400 

b 

VS+{0.03B)2 

15 

14 

13 

11 

10 

8 

7 


In the above formulae, 


BR{hi - 

invisible) 

= BR{hi ^ 

■ XiX\) + BR{h^ ■ 

aiai)[BR{ai 

(12) 

BR{h2 - 

invisible) 

= BR{h2^ 

• X?X?) 





-\-BR{h2 ■ 

—aiai)[i?-R(ai - 

- 




-\-BR{h2 ■ 

— hihi)[BR{hi - 

invisible)]'^ 

(13) 

BR{hs - 

invisible) 

= BR{hs^ 

• X?X?) 





-\-BR{h'j, ■ 

—aiai)[i?-R(ai - 

- 




-\-BR{h'^ ■ 

— hihi)[BR{hi - 

invisible)]'^ 




-\-BR{h^ 

hih2)BR{hi — 

^ invisible)BR{h 2 

— invisible) 



-\-BR{h'j, ■ 

—>■ h2h2)\BR{h2 - 

invisible)]'^. 

(14) 
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